#include "functions.h"

using namespace std;

int main()
{
    int n = 41;
    double** A = new double*[n];
    double** Q = new double*[n];
    for(int i = 0; i < n; i++)
    {
        A[i] = new double[n]();
        Q[i] = new double[n]();
    }
    for(int i = 1; i < n; i++)
    {
        A[i][i-1] = 1.0;
    }
    A[0][n-1] = -1.0;
    A[3][n-1] = -1.0;
    Implicit_QR(A, Q, pow(10, -5), n);
    return 0;
};